THE DISCRETE FOURIER TRANSFORM: 
A CANONICAL BASIS OF EIGENFUNCTIONS 
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ABSTRACT 

The discrete Fourier transform (DFT) is an important opera- 
tor which acts on the Hilbert space of complex valued func- 
tions on the ring Z/iVZ. In the case where N = p is an odd 
prime number, we exhibit a canonical basis <t> of eigenvectors 
for the DFT. The transition matrix © from the standard ba- 
sis to 4> defines a novel transform which we call the discrete 
oscillator transform (DOT for short). Finally, we describe a 
fast algorithm for computing © in certain cases. 

1. INTRODUCTION 

The discrete Fourier transform (DFT) is probably one of 
the most important operators in modern science. It is om- 
nipresent in various fields of discrete mathematics and en- 
gineering, including combinatorics, number theory, com- 
puter science and, last but probably not least, digital sig- 
nal processing. Formally, the DFT is a family {Fn} of 
unitary operators, where each Fn acts on the Hilbert space 
J#n= C (Z/iVZ) by the formula 

F N [f]{ w ) = ^= £ e%«f(t). 

Although, so widely used, the spectral properties of the 
DFT remains to some extent still mysterious. For example, 
the calculation of the multiplicities of its eigenvalues, which 
was first carried out by Gauss, is quite involved and requires 
a multiple of number theoretic manipulations OATl . 

A primary motivation for studying the eigenvectors of the 
DFT comes from digital signal processing. Here, a function 
is considered in two basic realizations: The time realization 
and the frequency realization. Each realization, yields in- 
formation on different attributes of the function. The DFT 
operator acts as a dictionary between these two realizations 

Time ^ Frequency. 

From this point of view, it is natural to look for a di- 
agonalization basis, namely, a basis of eigenvectors (eigen 
modes) for Fn- In this regard, the main conceptual difficulty 
comes from the fact that the diagonalization problem is ill- 
defined, since, Fn is an operator of order 4, i.e., FA = Id, 



which means that it has at most four eigenvalues ±l,±z, 
therefore each appears with large multiplicity (We assume 
JV>4). 

An interesting approach to the resolution of this diffi- 
culty, motivated from results in continuous Fourier analy- 
sis, was developed by Griinbaum in |G|. In that approach, a 
tridiagonal operator Sn which commutes with Fn and admits 
a simple spectrum is introduced. This enable him to give 
a basis of eigenfunctions for the DFT. Specifically, Sn ap- 
pears as a certain discrete analogue of the differential opera- 
tor D — d 2 — t 2 which commutes with the continuous Fourier 
transform. 

1.1 Main results of this paper 

In this paper we describe a representation theoretic approach 
to the diagonalization problem of the DFT in the case when 
N = p is an odd prime number. Our approach, puts to the 
forefront the Weil representation |PW1 of the finite symplectic 
group Sp = SLaffip) as the fundamental object underlying 
harmonic analysis in the finite setting. Specifically, we ex- 
hibit a canonical basis <P p of eigenvectors for the DFT. We 
also describe the transition matrix @ /7 from the standard basis 
to <t>p, which we call the discrete oscillator transform (DOT 
for short). In addition, in the case p = 1 (mod 4), we describe 
a fast algorithm for computing & p (FOT for short). 

It is our general feeling that the Weil representation yields 
a transparent explanation to many classical results in finite 
harmonic analysis. To justify this claim, we describe an alter- 
native method for calculating the multiplicities of the eigen- 
values for the DFT, a method we believe is more suggestive 
then the classical calculations. 

The rest of the introduction is devoted to a more detailed 
account of the main ideas and results of this paper. 

1.2 Symmetries of the DFT 

Let us fix an odd prime number p and for the rest of the intro- 
duction suppress the subscript p from all notations. Gener- 
ally, when a (diagonalizable) linear operator A has eigenval- 
ues admitting large multiplicities, it suggest that there exists 
a group G — GaC GL (J^) of "hidden" symmetries consist- 
ing of operators which commute with A. Alas, usually the 
problem of computing the group G is formidable and, in fact, 



equivalent to the problem of diagonalizing A. If the operator 
A arise "naturally", there is a chance that the group G can be 
effectively described. In preferred situations, G is commu- 
tative and large enough so that all degeneracies are resolved 
and the spaces of common eigenvectors with respect to G are 
one-dimensional. The basis of common eigenvectors with 
respect to G establishes a distinguish choice of eigenvectors 
for A. Philosophically, we can say that it is more correct to 
consider from start the group G instead of the single operator 
A. 

Interestingly, the DFT operator F = F p admits a natural 
group of symmetries Gf , which, in addition, can be effec- 
tively described using the Weil representation. For the sake 
of the introduction, it is enough to know that the Weil repre- 
sentation in this setting is a unitary representation p : Sp — > 
U (Jff) and the key observation is that F is proportional to a 
single operator p (w) . The group Gf is the image under p of 
the centralizer subgroup T w of w in Sp. 

1.3 The algebraic torus associated to the DFT 

The subgroup T w can be computed explicitly and is of a very 
"nice" type, it consists of rational points of a maximal al- 
gebraic torus in Sp which in plain language means that it 
is maximal commutative subgroup in Sp, consisting of ele- 
ments which are diagonalizable over some field extension. 
Restricting the Weil representation to the subgroup T w yields 
a collection Gf = {p (g) ■ g € T w } of commuting operators, 
each acts unitarily on the Hilbert space Jf? and commutes 
with F. This, in turn, yields a decomposition, stable under 
Fourier transform, into character spaces 

je=Qje x , (i) 

where % runs in the set of (complex valued) characters of 
r w , namely, if v G J#x then p (g) v = % (g) v. The main tech- 
nical statement of this paper, Theorem 3, roughly says that 
dim J4? x = 1 for every % which appears in £[). Choosing a 
unit representative (j) x £ ,yf x for every %, gives the canonical 
basis <t> = {0 Z } of eigenvectors for F. The oscillator trans- 
form © sends a function / e Jf? to the coefficients in the 
unique expansion 

The fine behavior of F and © is governed by the (split type) 
structure of T w , which changes depending on the value of 
the prime p modulo 4. This phenomena has several conse- 
quences. In particular, it gives a transparent explanation to 
the precise way the multiplicities of the eigenvalues of F de- 
pend on the prime p. Another, algorithmic, consequence is 
related to the existence of a fast algorithm for computing 0. 

1.4 Properties of eigenvectors 

The character vectors § x satisfy many interesting properties 
and are objects of study in their own right. A comprehensive 
treatment of this aspect of the theory appears in 1GHSI . 

1.5 Generalizations 

1.5.1 Field extensions 

All the results in this paper were stated for the basic finite 
field F p , for the reason of making the terminology more ac- 
cessible. In fact, all the results can be stated and proved for 



any field extension of the form F q , q = p". In most places 
one should only need to replace p by q. 

1.6 Structure of the paper 

We begin by discussing the finite Heisenberg group and the 
Heisenberg representation. Next we introduce the Weil rep- 
resentation of the finite symplectic group, first it is described 
in abstract terms and then more explicitly invoking the idea 
of invariant presentation of an operator. We proceed to dis- 
cuss the theory of tori in the one-dimensional Weil repre- 
sentation, we explain how to associate to a maximal torus 
T C SL2, a transform ©7- called the oscillator transform. We 
describe a fast algorithm for computing ©7- in the case T is a 
split torus. The theory is then applied to the specific torus as- 
sociated with the DFT operator. We finish with a treatment of 
the multiplicity problem for the DFT, from the representation 
theoretic perspective. 
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2. THE OSCILLATOR TRANSFORM 

2.1 The Heisenberg group 

Let (V,(0) be a two-dimensional symplectic vector space 
over the finite field F p . The reader should think of V as F p x 
F p with the standard form 00 ((t,w) , (Y,w')) = Tw' — wx'. 
Considering V as an abelian group, it admits a non-trivial 
central extension called the Heisenberg group. Concretely, 
the group H can be presented as the set H = V x F p with the 
multiplication given by 

(v,z) • (v',z!) = {v + v',z + z' + \(o(y,V)). 

The center of H is Z = Z(H) = {(0,z) : z G F p } . The 
symplectic group Sp = Sp(V,Co), which in this case is iso- 
morphic to SL2 (F p ), acts by automorphism of H through its 
action on the V-coordinate. 

2.2 The Heisenberg representation 

One of the most important attributes of the group H is that it 
admits, principally, a unique irreducible representation. The 
precise statement goes as follows. Let iff : Z — > C x be a char- 

acter of the center. For example we can take y(z) =e p . 
Theorem l.(Stone-von Neumann) 

There exists a unique (up to isomorphism) irreducible 
unitary representation {%, H,Ji?) with the center acting by 
i/A, i.e., 7Z\ Z = yz-Idj?. 

The representation % which appears in the above theorem 
will be called the Heisenberg representation. 



2.2.7 Standard realization of the Heisenberg representa- 
tion. 

The Heisenberg representation (7Z,H,Jf?) can be realized 
as follows: 3^ is the Hilbert space C(F p ) of complex val- 
ued functions on the finite line, with the standard Hermitian 
product. The action % is given by jt(t,0)[/] it) = f(t+ t), 
n(0,w)[f] (t) = yr{wt)f(?) and n{z)[f] (t) = y(z)f(t). We 
call this explicit realization the standard realization. 

2.3 The Weil representation 

A direct consequence of Theorem 1 is the existence of a 
projective representation p : Sp — > PGL(Jf). The construc- 
tion of p out of the Heisenberg representation n is due to 
Weil [W| and it goes as follows. Considering the Heisen- 
berg representation % and an element g 6 Sp, one can de- 
fine a new representation iz g acting on the same Hilbert 
space via n 8 (h) — % (g (h)). Clearly both % and n s have the 
same central character y/ hence by Theorem 1 they are iso- 
morphic. Since the space Hom//(7C, 7T A ') is one-dimensional, 
choosing for every g 6 Sp a non-zero representative p(g) £ 
Hom//(7T, 7Z g ) gives the required projective representation. In 
more concrete terms, the projective representation p is char- 
acterized by the formula 



p{g)%{h)p(g- 1 )=n( s {h)), 



(2) 



for every g £ Sp and h £ H. It is a peculiar phenomenon of 
the finite field setting that the projective representation p can 
be linearized into an honest representation. 
Theorem 2 There exists a uniqu$\ linear representation 

p:Sp — >GL{3f), 

which satisfies equation (0. 

2.3.1 Invariant presentation of the Weil representation 

Let us denote by C (H, \\f) the space of (complex valued) 
functions on H which are y/-equivariant with respect to the 
action of the center, namely, a function / £ C (H, y/) satis- 
fies / (zh) = yf(z)f (h) for every z £ Z, h £ H. Given an 
operator A £ End (Jtf), it can be written in a unique way as 
A = % (Ka), where K A £ C (H, ty/^ 1 ) and % denotes the ex- 
tended action n (Ka) — L Ka (h) it (h) . The function Ka is 
hen 

called the kernel of A and it is given by the matrix coefficient 



K A (h) 



1 



dim Jf 



7V (Arc (/T 1 )) 



(3) 



In the context of the Heisenberg representation, formula 
(01 is usually referred to as the Weyl transform. Using the 
Weyl transform one is able to give an explicit description 
of the Weil representation. The idea 1GH1II is to write 
each operator p (g), g £ Sp in terms of its kernel function 
Kg = Kpig) £ C (H, V -1 )- The following formula is taken 
from [GH1 1 



K g (v,z) 



-1) 



dimJf 



a (det (k (g) +/)) y (I (o (k (g) v, v) + z) 

(4) 



for every g £ Sp such that g — I is invertible, where a de- 
notes the unique quadratic character (Legendre character) of 
the multiplicative group and k is the Cayley transform 

2.4 The theory of tori 

A maximal (algebraic) torus in Sp is a maximal commutative 
subgroup which becomes diagonalizable over some field ex- 
tension. There exists two conjugacy classes of maximal (al- 
gebraic) tori in Sp. The first class consists of those tori which 
are diagonalizable already over ¥ p or equivalently those are 
the tori that are conjugated to the standard diagonal torus 



A = 











:a£¥ p 



A torus in this class is called a split torus. The second class 
consists of those tori which become diagonalizable over a 
quadratic extension ¥ p 2 or equivalently those are tori which 
are not conjugated to A. A torus in this class is called a non- 
split torus (sometimes it is called inert torus). 
Example l(Example of a non-split torus) 

It might be suggestive to explain further the notion of 
non-split torus by exploring, first, the analogue notion in 
the more familiar setting of the field R. Here, the stan- 
dard example of a maximal non-split torus is the circle group 
SO{2) C SZ^R)- Indeed, it is a maximal commutative sub- 
group which becomes diagonalizable when considered over 
the extension field C of complex numbers. The above anal- 
ogy suggests a way to construct an example of a maximal 
non-split torus in the finite field setting as well. 

Let us identify the symplectic plane V = ¥ p x ¥ p with 
the quadratic extension ¥ p 2. Under this identification, ¥ 2 
acts on V and for every g £ ¥ p i we have 0)(gu,gv) — 
det (g) CO (u,v), which implies that the group 



Tns ={#eF p x 2 : det (g) = l} 



naturally lies in Sp. The group T ns is an example of a non- 
split torus which the reader might think of as the "finite cir- 
cle ". 

2.4.1 Decompositions with respect to a maximal torus 

Restricting the Weil representation to a maximal torus T C 
Sp yields a decomposition 



x 



(5) 



where % runs in the set T v of complex valued characters 
of the torus T . More concretely, choosing a generatofl 
t £ T, the decomposition (0 naturally corresponds to the 
eigenspaces of the linear operator p (t). The decomposition 
(0 depends on the split type of T. Let aj denote the unique 
quadratic character of T. 
Theorem 3 (|GH2|) IfT is a split torus, then 



X = Ot- 



1 Unique, except in the case the finite field is F3 . For the canonical choice 
in the latter case see f GHll . 



2 A maximal torus T in SL,2(¥ p ) is a cyclic group, thus there exists a 
generator. 



If T is a non-split torus, then 

x \ X = or- 

2.5 The discrete oscillator transform associated to a 
maximal torus 

Let us fix a maximal torus T . Every vector v G j^ 5 can be 
written uniquely as a direct sum v = £v x with G and 
% runs in / = Spec r (J$?) - the spectral support of Jf? with 
respect to T consisting of all characters % G T v such that 
dim J4? x 0. Let us choose, in addition, a collection of unit 
vectors <j> x G % G 7 and let = Lfe- We define the 
transform ® r = : C(7) by @r [v] (#) = (v,^). 

We will call the transform ©7- the discrete oscillator trans- 
form (DOT for short) with respect to the torus T and the test 
vector . 

Remark 1 We note that in the case T is a non-split torus, 
©r maps ffl isomorphically to C(I). In the case T is a 
split torus, ©r has a kernel consisting of v G Jf such that 

(v,<M=o. 

2. 5. 7 77ie oscillator transform ( integral form ) 

Let J( r : C(T) -> C(r v ) denote the Mellin transform 

1 



■^r [/](*) 



for / G C (r). Let us denote by otj : — * C (T) the matrix 
coefficient mj [v] (g) = (v,p (g -1 ) 0) for v G J^. 
Lemma 1 (EH2)) We /zave 



7 



lj o my. 



2.5.2 Twf oscillator transforms 

In practice, it is desirable to have a "fast" algorithm for com- 
puting the oscillator transform (FOT for short). We work 
in the following setting. The vector v is considered in the 
standard realization M 1 = C (F p ) (see 12.2. it , in this context 
the oscillator transform gives the transition matrix between 
the basis of delta functions and the basis {^>x} °f charac- 
ter vectors. We will show that when T is a split torus and 
for an appropriate choice of 0, the oscillator transform can 
be computed in O (plog(p)) arithmetic operations. Princi- 
pally, what we will show is that the computation reduces 
to an application of DFT followed by an application of the 
standard Mellin transform, both transforms admit a fast al- 
gorithm JCT) . 

Assume T is a split torus. Since all split tori are con- 
jugated to one another, there exists, in particular, an ele- 
ment s G Sp conjugating T with the standard diagonal torus 
A. In more details, we have a homomorphism of groups 
Ad s -T — * A sending g G T to Ad s (g) = sgs' 1 G A. Dually, 
we have a homomorphism Ad^ :A V ~^T V between the cor- 
responding groups of characters. 

The main idea is to relate the oscillator transform with 
respect to T with the oscillator transform with respect to A. 
The relation is specified in the following simple lemma. 
Lemma 2 ( llGH2l ) We have 



Remark 2 Roughly speaking, (0 means that (up to a 
"reparametrization" of T v by A v using Ac/ 4 V ) the oscillator 
transform of a vector v G with respect to the torus T is 
the same as the oscillator transform of the vector p (s) v with 
respect to the diagonal torus A. 

In order to finish the construction we need to specify two 
basic facts about the Weil representation in the standard re- 
alization. 

• The standard torus A acts by (normalized) scalings, the 

f a 



precise formula of p (g) for g = I q j G A is 

P{g)\f] (*) = ff(a)/(<w). 

• Every operator p (g) , g G Sp, can be written in the form 
p (g) = M g[ oF o M g o S a where M gl ,M g . are the opera- 
tors of multiplication by some functions g\,g2 G C (F p ), 
S a is the operator of scaling by a G F* and F is the DFT. 

p\f}(y) = J 7= L v(yx)f(x). 

VP xe¥ p 

Given a function / G C (F p ), applying formula © with 
<p = p (s)^ 1 8i yields 

®^[/](A £ / J v (z)) = ^- £ a(a)x(a)p(s)[f](a), 

P «GF* 

(7) 

for every % G A v . In conclusion, formula © implies that 
®T,d> [f] can be computed by, first, applying the operator p (s) 
to / and then applying Mellin transform to the result. 
Problem 1 

Does there exists a fast algorithm for computing the os- 
cillator transform associated to a non-split torus? 

2.6 Diagonalization of the discrete Fourier transform 

In this subsection we apply the previous development in or- 
der to exhibit a canonical basis of eigenvectors for the DFT. 
We will show that the DFT can be naturally identified (up 
to a normalization scalar) with an operator p (w) in the Weil 
representation, where w is an element in a maximal torus 
7w C Sp. We take wG Sp = SLq (F p ) to be the Weyl element 



w = 



1 

-1 



Lemma 3([GH2]) We have 



F = C-p(w), 



.2-1 



(Aj/)*o© r ,«,=© 



A,p(s)<j> 



>p(j). 



(6) 



where C = i 

Lemma 3 implies that the diagonalization problems of the 
operators F and p (w) are equivalent. The second problem 
can be approached using representation theory, which is what 
we are going to do next. 

Let us denote by T w the centralizer of w in Sp, namely T w 
consists of all elements g G Sp such that gw=wg, in particu- 
lar we have that wG T w . 

Proposition 1 ([GH2)) The group T w is a maximal torus. 
Moreover the split type of T w depends on the prime p in the 



following way: T w is a split torus when p = 1 [mod 4) and is 
a non-split torus when p = 3 (mod 4). 

Proposition 1 has several consequences. First conse- 
quence is that choosing a unit character vector <p x E Jt? x for 
every % £ Spec rw (Jtf) gives a canonical (up to normaliz- 
ing unitary constants) choice of eigenvectors for the DFT0. 
Second, more mysterious consequence is that although the 
formula of the DFT is uniform in p, its qualitative behav- 
ior changes dramatically between the cases when p = 1 (mod 
4) and p = 3 (mod 4). This is manifested in the structure 
of the group of symmetries: In the first case, the group of 
symmetries is a split torus consisting of p — 1 elements and 
in the second case it is a non-split torus consisting of p + 1 
elements. It also seems that the structure of the symmetry 
group is important from the algorithmic perspective, in the 
case p = 1 (mod 4) we built a fast algorithm for computing 
©, while in the case p = 3 (mod 4) the existence of such an 
algorithm remains open (see Problem 1). 

2.6.1 Multiplicities of eigenvalues of the DFT 

Considering the group T w we can give a transparent compu- 
tation of the eigenvalues multiplicities for the operator p (w). 
First we note that, since w is an element of order 4, the eigen- 
values of p (w) lies in the set {±1,±/}. For A S {±1,±/}, let 
mi denote the multiplicity of the eigenvalue A . We observe 
that 

m X = dim J^, 

where consists of all characters % <= Spec rw (Jf?) such that 
^(w) = A. The result now follows easily from Theorem 3, 
applied to the torus T w . We treat separately the split an non- 
split cases. 

• Assume T w is a split torus, which happens when p = 
1 (mod 4), namely, p = 41 + 1, I G N. Since dim^. = 1 

for x 7^ Or w it follows that m ±i = ^ = I. We are left to 
determine the values of m±\, which depend on whether 

°Zw (w) =wT is 1 or —1. 

Since w is an element of order 4 in T w we get that 



Summarizing, the multiplicities of the operator p (w) are 



<*r w (w) = 



1, p = l(mod8), 
-1, p = 5(mod8), 



which implies that when p = 1 (mod 8) then m\ =1 + 1 
and m_i = / and when p = 5(mod 8) then m\ = / and 

OT_ l =1+1. 

Assume T w is a non-split torus, which happens when p = 
3 (mod 4), namely, p = 41 + 3, I e N. Since dim^ = 1 

for x 7^ Or w it follows that m±, = = / + 1 . The values 

of m±i depend on whether &t w (w) =wT~ is 1 or —1. 
Since w is an element of order 4 in T w we get that 





mi 


m_i 


mj 


m-i 


p = 8k+l 


2k+\ 


2k 


2k 


2k 


p = 8k + 3 


2k 


2k+l 


2k+l 


2k+\ 


p = U + 5 


2k+l 


2k + 2 


2k+l 


2k+l 


p = $k + l 


2k + 2 


2k+l 


2k + 2 


2k + 2 



(8) 



Considering now the DFT operator F. If we denote by 
rijx, fl € {±1,±/} the multiplicity of the eigenvalue jj. of F 
then the values of can be deduced from table[8]by invok- 

p-i 

ing the relation n^ = mx where A = i^~ ■ pL (see Lemma 3). 
Summarizing, the multiplicities of the DFT are 







n-\ 


Hi 


n-i 


p 


= 41 + 1 


1 + 1 


I 


I 


I 


p 


= 41 + 3 


l + l 


l + l 


l + l 


I 



For a comprehensive treatment of the multiplicity problem 
from a more classical point of view see II ATI . Other applica- 
tions of Lemma 3 appear in BGHH1 . 
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°T V (w) = 



1, /? = 7(mod8), 
-1, /? = 3(mod8), 



which implies that when p = 7 (mod 8) then m\ = I and 
m_ i = / + 1 and when p = 3 (mod 8) then m\ =1 + 1 and 
m_i = I. 



3 In the case T w is a split torus there is a slight ambiguity in the choice of 
a character vector with respect to <Tp . This ambiguity can be resolved by 
further investigation which we will not discuss here. 



